A Study of Synthetic and Observed Ha Spectra of TT Hydrae 
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ABSTRACT 

The formation and properties of accretion discs and circumstellar material in Algol-type 
systems is not very well understood. In order to study the underlying physics of these structures, 
we have calculated synthetic Ha spectra of TT Hya, which is an Algol-type eclipsing binary with 
an accretion disc. Both the primary and secondary stars were considered in the calculations as 
well as a disc surrounding the primary. The Roche model for the secondary star was assumed. 
The synthetic spectra cover all the phases including primary eclipse and are compared with the 
observed spectra. The influence of various effects and free parameters of the disc on the emerging 
spectrum was studied. This enabled us to put some constraints on the geometry, temperature, 
density and velocity fields within the disc. Differences found between the observed and synthetic 
spectra unravel the existence of a gas stream as well as a hotter disc-gas interaction region. 
An additional cooler circumstellar region between the CI and C2 Roche surfaces is suggested 
to account for various observed effects. A new computer code called SHELLSPEC was created 
for this purpose and is briefly described in this work as well. It is designed to solve simple 
radiative transfer along the line of sight in 3D moving media. The scattered light from a central 
object is taken into account assuming an optically thin environment. Output intensities are then 
integrated through the 2D projection surface of a 3D object. The assumptions of the code include 
LTE and optional known state quantities and velocity fields in 3D. 
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Stars: 



TT Hydrae (HD 97528, HIP 54807, SAO 
179648, V = 7.27™, a = = -26°28') is 

a typical Algol-type eclipsing binary system with 
an orbital period of Porb — 6.95 days (Kulkarni & 
Abhyankar 1980) and consists of a hotter B9.5 V 
main sequence primary and a cooler evolved Kl 
III-IV secondary star filling its Roche lobe. There 
is evidence of the circumstellar material in the 
form of an accretion disc (Plavec & Polidan 1976) 
surrounding the primary as well as an indication 
of a gas stream (Peters & Polidan 1998). Kulkarni 
& Abhyankar (1980) obtained UBV photometric 
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observations and photometric elements of the sys- 
tem. These observations were re-examined by 
Etzel (1988) who found from photometry, spec- 
trophotometry and spectroscopy: Toff — 9800K, 
wsini = 168 ± 5 kms^^ for the primary, and 
Tcff = 4670 - 4850K for the secondary. He also 
suggested, based on the Inglis-Teller formula, that 
the gas producing the UV excess and Balmer line 
emission has an electron number density lower 
than 10^^ cm"'^. He also estimated a distance of 
193 pc to the binary. Eaton & Henry (1992) esti- 
mated wsini = 43 ± 3 kms~^ for the secondary. 
Plavec (1988) studied the lUE spectra of TT Hya 
and concluded from the presence of Fe n absorp- 
tion lines that the vertical dimension of the accre- 
tion disc is significant and is at least comparable 
to the diameter of the primary. He estimated a 
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distance to the system of 194 pc, and pointed out 
that the presence of super- ionized emission hnes of 
Si IV, C IV, N V poses a question about the source 
of ionization. Peters (1989) estimated from the 
echpses of the Ha emission region that a fairly 
symmetrical disc fills up to 95% of the Roche lobe 
radius. Moreover, she argued that the disc must 
be rather flat with the vertical dimension com- 
parable or less than the diameter of the primary. 
The mean electron number density in the disc was 
estimated to be ~ 10^° cm~'^. The depth of the 
Ha absorption was found to be strongly variable 
with phase and was deepest shortly before and 
after primary eclipse; a feature also observed in 
other Algols (e.g., Richards 1993). Peters (1989) 
also suggested that the additional absorption in 
the center of the Ha line probably comes from the 
dense inner region of the accretion disc. She also 
observed that the Ha core is blueshifted relative 
to the primary by about 50 km s~^ and interpreted 
it as evidence of the mass outflow. 

Sahade & Cesco (1946) determined the orbital 
parameters of the primary and identified an ec- 
centric orbit from the CaiiK and H lines. Pop- 
per (1989) measured radial velocities of the pri- 
mary and secondary and found that the motion 
of the secondary was consistent with the circular 
orbit {K = 132 kms~^) but not in the case of 
the primary. Van Hamme & Wilson (1993) re- 
analyzed the light and velocity curves separately 
and simultaneously and used a physical model 
to obtain the following parameters: separation 
a = 22.63 ± O.12i?0, effective temperature of the 
primary Ti = 9800K, secondary Ta = 4361K, 
radius of the primary i?i = 1.95i?Q, secondary 
i?2 = 5.87i?0, mass of the primary Mi = 2.63M0, 
secondary M2 ~ O.59A'/0, and mass ratio q = 
0.2261 ± 0.0008. Albright & Richards (1996) ob- 
tained the first Doppler tomogram to reveal an 
indirect image of the accretion disc. Hipparcos 
(ESA 1997) obtained a parallax of 6.50±0.95 mas. 
Peters & Polidan (1998) observed red-shifted ab- 
sorption in Ni, Nil lines in their FUV spectra 
at the phase 0.95 and interpreted it as evidence 
of the gas stream and inferred the rate of inflow 
to be greater than lO~"'^^M0yr~-'^. Most recently, 
Richards & Albright (1999) studied the Ha differ- 
ence profiles and properties of the accretion region 
and put TT Hya into the context of other Algol 
binaries. 



It would be interesting to verify the above es- 
timates, conclusions, and interpretations on more 
physical grounds using some more sophisticated 
methods involving synthetic spectra which would 
take into account the circumstellar matter. Un- 
fortunately, there are not many tools available for 
such a task at present. 

There are sophisticated computer codes for cal- 
culating and inverting light curves or spectra of bi- 
nary stars with various shapes or geometry includ- 
ing the Roche model (Lucy 1968; Wilson & Devin- 
ney 1971; Mochnacki & Doughty 1972; Rucin- 
ski 1973; Hill 1979; Zhang et al. 1986; Djura- 
sevic 1992; Drechsel et al. 1994; Vinko et al. 
1996; Hadrava 1997; Bradstreet fc Steelman 2002; 
PribuUa 2004). In these codes, the stars are as- 
sumed to be nontransparent and stripped of any 
circumstellar matter. An interesting approach to 
study extended semitransparent atmospheres of 
some eclipsing binary components was developed 
by Cherepashchuk et al. (1984). Often the 3D 
model of the circumstellar matter (behavior of 
state quantities and velocity field) is known or ex- 
pected as a result of hydrodynamic simulations or 
observational constraints (see for example Karet- 
nikovctal. 1995; Richards & Ratliff 1998). Unfor- 
tunately, 3D Non-LTE (NLTE) calculations and 
spectral synthesis including complex hydrodynam- 
ics are difficult to carry out so one alternative has 
been to perform a simple volume integration of 
emissivity, which is often too oversimplified for the 
particular problem. 

On the other hand, highly sophisticated model 
atmospheres and spectrum synthesis codes have 
been developed assuming NLTE and plane-parallel 
atmospheres of hot stars (Hubeny 1988; Hubeny & 
Lanz 1992, 1995; Hubeny et al. 1994), spherically 
symmetric atmospheres (Kubat 2001), or stellar 
winds (Krticka & Kubat 2002). There are also 
sophisticated stationary plane-parallel line-blan- 
keted model atmospheres and spectrum synthesis 
codes for a large variety of stars assuming LTE 
(Kurucz 1993a, 1993b; Smith & Dworetsky 1988; 
Piskunov 1992; and many others). However, these 
are very specialized codes and their main purpose 
is to calculate the spectrum emerging from a stel- 
lar atmosphere, and it is difficult to apply them 
to the various cases outlined above. An excep- 
tion is the special case of circumstellar matter in 
the form of accretion discs in cataclysmic variables 
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(CVs). In this case, the disc is either approxi- 
mated by a set of geometrically thin, but optically 
thick, static local atmospheres and the output ra- 
diation is a sum of properly Doppler shifted local 
emerging intensities (Orosz & Wade 2003; Wade & 
Hubeny 1998; la Dous 1989) or, in case of optically 
thin discs or accretion disc winds, the Sobolev ap- 
proximation is used (Proga et al. 2002; Long & 
Knigge 2002; Rybicki & Hummer 1983). Linnell 
& Hubeny (1996) developed a package of codes 
which can calculate light curves or spectra of inter- 
acting binary stars including optically thick (non- 
transparent) disc. However, it is often necessary 
to solve the radiative transfer in a moving disc at 
least along the line of sight as demonstrated by 
Horne & Marsh (1986). Nevertheless, their calcu- 
lations assumed a linear shear and neglected stim- 
ulated emission, continuum opacity, scattering, as 
well as a central star. 

In this paper, we attempt to bridge the gap 
in the previously mentioned approaches using our 
new code SHELLSPEC. It is a tool which solves in 
LTE the simple radiative transfer along the line 
of sight in an optional optically thin 3D moving 
medium. Optional transparent (or non transpar- 
ent) objects such as a spot, disc, stream, jet, shell 
or stars as well as an empty space may be defined 
(or embedded) in 3D and their composite syn- 
thetic spectrum calculated. The stars may have 
the Roche geometry and known intrinsic spectra. 
The code is quite a multi-purpose, independent, 
and flexible tool which can also calculate a light 
curve or a trailing spectrogram where necessary, 
or can be used to study various objects or effects. 
In Section 2, we describe the basic astrophysics in- 
cluded in the code. In Section 3, we briefly men- 
tion numerical methods used in the code. In Sec- 
tion 4.1, we study in more detail the effects of var- 
ious input parameters on the spectrum of a star 
which has a circumstcllar disc. Budaj & Richards 
(2004) made preliminary calculations for a single 
star with a disc with parameters similar to TT 
Hya. However, these were based on the stellar 
parameters found by Etzel (1988), a wedge shape 
geometry of the disc with a relatively large inner 
radius (2.5i?0), and no limb darkening of the pri- 
mary star. More realistic stellar parameters by 
Van Hamme & Wilson (1993) and a slab shaped 
disc geometry will be used in this work. Then in 
Section 4.2, we will make more detailed calcula- 



tions of the system for different phases involving 
both stars and a disc and will compare our new 
synthetic spectra to real observations of TT Hya. 

2. Basic Astrophysics 

This section describes the basic astrophysics in- 
cluded in the code. For more detailed information 
and a manual for the code, the reader is referred 
to Budaj k Richards (2004). 

2.1. Radiative Transfer 

In the following analysis, the calculations are 
carried out in the observer's Cartesian frame with 
z pointing towards the observer. The radiative 
transfer equation along the line of sight is: 



(1) 



where is the specific monochromatic intensity 
at the frequency v, Xv is the opacity, e,j is the 
emissivity and z is the distance along the beam. 
It is convenient to split the opacity into two contri- 
butions, the true absorption and the scattering 

Xu = l^v + CFv (2) 

Assuming LTE, the line opacity corrected for a 
stimulated emission is stated simply as: 

^Irne ^ (1 _ g"^ )7V;B,„/liy^,„ (l. - zyo)(47r)-l (3) 

where h is the Planck constant, hv is the energy 
of the transition from the lower level I to the up- 
per level u, fc is the Boltzmann constant, T is the 
temperature, A^; is the population of the ?-th state 
of the corresponding ion, and Biu is the Einstein 
coefficient for the whole solid angle. The velocity 
field enters the equation via the shifted normalized 
Voigt profile y;„(i' — vq) where 



v^iz) 



(4) 



where viu is the laboratory frequency of the line 
and Vz = v.n is the radial velocity (positive to- 
wards the observer) or projection of the local ve- 
locity vector V to the line of sight unit vector n. 
The Einstein coefficient, Biu, is related to the os- 
cillator strength, /;„, by: 



Biu 



mechviu 



(5) 
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where e, mo are the electron charge and mass, re- 
spectively and c is the speed of light. The shape of 
the Voigt profile is determined by the thermal and 
the microturbulent broadening, ^^^6) characterized 
by the Doppler half-width 



Al/D = - 



2kT 



+ v. 



trb 



(6) 



as well as by the damping broadening character- 
ized by the frame damping parameter 

a = 7/(47rAj/£,) (7) 
where the damping constant 

7 = iNat. + IStark + IVDW (8) 

includes the contribution from the Natural, Stark 
and Van der Waals broadening. In the case of 
LTE, all the line opacity is due to the true ab- 
sorption process i.e.: 



(9) 



We also included four other continuum opac- 
ity sources: the HI bound-free opacity, the 
HI free-free opacity, Thomson scattering and 
Rayleigh scattering on neutral hydrogen, denoted 
by nff^''^, nff^^^, cr"^^, a^^ , respectively. 
The total true absorption k^, is the sum of the 
three opacity sources: 



^ _ Jine , HIbf , Hlff 



The total scattering is: 



^ — ^TS j_ RS 



(10) 



(11) 



The thermal emissivity associated with the true 
absorption can then be written as 



(12) 



where is the Planck function for the local value 
of the temperature. For scattering emissivity we 
have 

= JJ <j{v',v,n\n)I{v',n')di,''^ (13) 

where (j{v' ,v,n' ,n) is the scattering coefficient 
containing the general redistribution function. It 
is this term which causes the main difficulty, since 



apart from redistributing the frequencies [v' v), 
it also couples the radiation in one direction n 
with the radiation field in all other directions n' . 
However, in many applications (e.g., optically thin 
shells) this term can either be neglected or gov- 
erned by the scattering of light from the central 
object. Wc assume coherent isotropic scattering 
(as seen from the scattering particle frame) from 
a blackbody or from a central spherical star with 
prccalculated surface intensity /* or flux F*. In 
this case the emissivity reduces to: 



SC T 



(14) 



where is the mean intensity. Ignoring limb 
darkening, J^, can be approximated by: 



(15) 



where lo is the solid angle subtended by the central 

star and 



a;/47r = ^ I 1 



r 



(16) 



where i?* is the radius of the central star and r is 
the distance from the center (of the star/grid) and 



1^1 



and 



Vi 



r.{v - Vi,) 



(17) 



(18) 



where v is the velocity field vector at the given 

point specified by the vector r and i;* is the ve- 
locity of the center of mass of the central object. 
For i?*/r << 1, an approximation including the 
limb darkening and the non-isotropic dipole phase 
function g{n',n) = |(1 + (n'.n)^) can be used: 



^(^l + ^la.J. 



where 



where 



F* 



p2 



47r 



F* =^1* 



(-1) 



where u is the limb darkening coeSicient. 
The total emissivity is then 



(19) 

(20) 
(21) 

(22) 
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and the total source function is: 



(23) 



The radiative transfer equation along the line 
of sight is solved (see Section 3) and the flux, F^, 
from the object at the Earth is then obtained by 
the integration of the output intensities through 
the 2D projection surface of the 3D object. 

2.2. Roche Geometry 

Both objects, star and companion, may have 
shapes according to the Roche model for detached 
or contact systems. Descriptions of the Roche 
model can be found in Kopal (1959), Limber 
(1963), Plavec & Kratochvil (1964), Mochnacki & 
Doughty (1972), Hilditch (2001) and many other 
papers and books. Most of them use cylindrical 
coordinates but wc will assume a Cartesian coor- 
dinate system (x, y, z) centered on one of the stars 
(labeled as 1) such that the companion (labeled as 
2) is at (1,0,0) and revolves around the z axis in the 
direction of the positive y axis. Let the mass ratio, 
q, always be m2/mi or 'companion/star' and g < 1 
will indicate the central star is heavier while > 1 
means the central star is lighter. We will also as- 
sume synchronous rotation i.e., the stellar surface 
rotates with the orbital period. Then, the normal- 
ized Roche potential, C(a;, y, z), is expressed as: 



C 



2q 



(l-hg)ri (l + g)r2 



where 



and 



ri = \/x^ +y'^ + z^, 
ri = V(a; - 1)2 + 2/2 + 



(24) 
(25) 

(26) 

The Roche surface of a detached component is de- 
fined as an cquipotcntial surface Cg = C{xa, yg. z^) 
passing through the sub-stellar point {xs,ys,Zs) 
(point on the surface of the star in between the 
stars, < Xs < l,j/s ~ Zs — 0) which is localized 
by the 'fill-in' parameter fi < 1. We define this 
by: 

fi = x,/L,,, fi = {l-Xs)/il-L,,) (27) 

for the primary and the secondary, respectively. 
Lix is the X coordinate of the LI point at 
(Lix,0, 0). The Roche equipotential surface Cg 



of a contact system will be defined by the fill-out 
parameter 1 < /o < 2: 



fo 



Cl-Cs 
CI-C2 



+ 1 



(28) 



where potentials C1,C2 correspond to the po- 
tentials at the LI and L2 points, respectively. 
First, we calculate LI, L2, Cg and x-boimdaries 
of the object using the Newton-Raphson iteration 
method in the x direction, and then the 3D shape 
of the surface is solved using the Newton-Raphson 
iteration in y and z coordinates with the precision 
of about 10"^. 

Gravity darkening is taken into account by 
varying the surface temperature according to the 
following law (Claret 1998): 



T/Tp = (g/gp) 



/3 



(29) 



where g is the normalized surface gravity, /3 is the 
gravity darkening exponent, Tp, and gp are the 
temperature and gravity at the rotation pole. The 
normalized gravity vector is g = {Cx,Cy, Cz) and 



where: 



g = ^Cl + Cl + Cl 



5C _dC ^ _dC 
OX ay oz 



(30) 



(31) 



The gravity darkening factor of the surface inten- 
sity is then calculated as: 



foD = B,{T)/B,{Tp). 



(32) 



Notice, that there is an imminent singularity in the 
calculations in the vicinity of LI, L2 points since 
gravity falls to zero which drags temperatures (a 
denominator in many equations) to zero. We avoid 
the problem by setting the lowest possible value 
of g/9p = 10^^. Limb darkening is taken into 
account using Eq. 42 and by calculating the cosine 
of the angle between the line of sight unit vector 
^ = {f^x^nyjiiz) and a normal to the surface: 



cos^ 



-n.g/g = , ■ (33) 



The reflection effect is not included in the present 
version. 
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3. Numerical Performance 

3.1. Solving the Radiative Transfer Equa- 
tion 

A number of optional objects (including trans- 
parent, nontransparent, empty space) can be de- 
fined within the model and the line of sight may 
cross more of them within a few grid points. A 
simple and stable method is needed to cope with 
such velocity, density, temperature fields which are 
optional and are allowed to be non-continuous. 
Our problem is that of the integration of a first 
order ordinary differential equation with known 
initial values at one boundary. The equation of 
the radiative transfer along the line of sight at the 
frequency v is integrated with an analogue of the 
mid-point method (Budaj & Richards 2004). It 
can be written in the discretised form: 



h+i 



Zi+l — Zj 



= ei+1/2 - Xi+i/2{Ii+i + (34) 



and simply integrated via the following recurrent 
formula: 

_ ^^+1/2 A -1/2 



where 



and 



where 



Xi+l/2iZi+l - Zi) 



(36) 



i+1/2 



'^»+l/2 
Xi-l-1/2 



, h=0 (37) 



Xi+1/2 - ^ • (3»j 



If the line of sight happens to hit a nontrans- 
parent object, the integration starts on the other 
side of the object with the boundary condition: 



or 



h = I*{v2)fLDfGD 



h = Bv^{Tes)fLDfGD 



(39) 



(40) 



where V2 is the Doppler shifted frequency in the 
comoving frame corresponding to the frequency v 
of the observer's frame: 



(41) 



and V* is the radial velocity of the surface of the 
nontransparent object where it intersects the line 
of sight. I* is the surface intensity of the non- 
transparent object perpendicular to the surface in 
the comoving (frozen to the surface) frame. If /* is 
not available it can be calculated e.g., using Eq. 21 
from the known surface flux. Rotation of the non- 
transparent objects is fully taken into account here 
by including it into the calculations of v*. /ld is 
a limb darkening factor: 



fLD = l — u + u cos 9 



(42) 



where u is the limb darkening coefficient and 9 is 
the angle between the normal to the surface of the 
nontransparent object and the line of sight. fcD 
is a gravity darkening factor which is important in 
the case of Roche geometry (see the Section 2.2) 
and in which case T^ff = Tp is the temperature at 
the rotation pole of a detached star or the tem- 
perature at the rotation pole of the more massive 
star in the case of a contact system, otherwise, it is 
the common effective temperature of the spherical 
star. If the line of sight happens to pass through 
an empty space this region is skipped and the in- 
tegration continues with I^+i = ij. 

The SHELLSPEC code enables the user to look 
on the object from different points of view and 
to calculate the corresponding spectra. The input 
model of the shell is defined in its 'body frozen' 
Cartesian coordinates {x" , y", z") with the z" axis 
corresponding to the intrinsic rotation axis of the 
model. The spectrum is always calculated in the 
observer's 'line of sight' Cartesian frame {x,y,z) 
with z pointing to the observer and which has the 
same center of coordinates. 

3.2. Additional Information and Adopted 
Routines 

The code assumes the known behavior of state 
quantities: temperature T, density p, and elec- 
tron number density ng. Solar abundances are 
assumed (Grevesse k. Sauval 1998) but the user 
is allowed to change all the element abundances. 
The level populations are obtained from the Boltz- 
mann and Saha equations. Partition functions 
(routine 'pfdwor') were taken from the UCLSYN 
code (Smith & Dworetsky 1988; Smith 1992). A 
FORTRAN 7 7 code containing the partition func- 
tion routines is also available in Budaj, Dworetsky 
& Smalley (2002). The Gaunt factors (routines 
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'gaunt', 'gfree') and Voigt profile (routine 'voigtO') 
are calculated with the subroutines taken from the 
SYNSPEC code (Hubeny et al. 1994). If the damp- 
ing constants arc not known they are estimated 
in the code in the way analogous to the SYNSPEC 
code. We also adopted the atomic data for chem- 
ical elements (routine 'stateO') from SYNSPEC and 
routines 'locate', 'hunt' from Numerical Recipes 
(Press et al. 1986). We also used a few sections 
from our previous original codes for calculations 
of radiative accelerations in stellar atmospheres of 
hot stars (Budaj & Dworetsky 2002). Apart from 
the above, the code was written from scratch and 
provides an independent tool to study a large va- 
riety of objects and effects. CGS units are used 
within the code and the manuscript, if not speci- 
fied otherwise. 

4. Application to Accretion Disc in TT 

Hya Type Systems 

4.1. The Effects of Various EVee Parame- 
ters on the Emerging Spectrum 

Before we modeled the complex spectrum of the 
two stars and the disc, we studied the effects of 
various free parameters on the spectrum. To com- 
prehend the problem more easily we restricted our 
calculations to the following assumptions. We ne- 
glected the secondary and studied only the pri- 
mary star surrounded by a disc. The primary 
was treated as a spherical blackbody with the fol- 
lowing parameters: effective temperature Toff = 
9800K, mass M = 2.63M0, radius R = IMRq 
(Van Hamme & Wilson 1993) and limb darken- 
ing u = 0.5. The disc was characterized by a 
slab of vertical half-width a = 2Rq which was 
further constrained by two spherical surfaces with 
the inner radius = 2Rq and the outer radius 
Rout = lOi?0. The electron number density of 
the disc was set to rig = 3 x 10*^ cm^'^, density 
p = 5 X IQ-is gcm-3, temperature T = 7000K, 
microturbulence Vtrb = kms~^, and inclination 
i = 82.84°. We adopted the following velocity field 
V of the disc: 



u}{r) = Y G— V = uj X r (43) 

where M is the mass of the central object, G is the 
gravitational constant, u>,r are the angular veloc- 
ity and radius vectors, respectively. The above 



prescription for the velocity reduces to the well 
known Keplerian velocity field in the limit of zero 
vertical thickness of the disc. The calculations 
were performed on the Cartesian 3D grid with the 
typical resolution in the disc plane of 121 x 121 
points. Figures 1 to 3 illustrate the effects of vary- 
ing just one of these free parameters of the system. 

The top of Figure 1 displays the effect of in- 
creasing the vertical half thickness of the disc from 
0.5 to 5i?0. The overall emission increases as we 
increase the emitting space volume, and the cen- 
tral absorption decreases mainly for the thicker 
discs. This occurs because the central depression 
is partly produced by the material projecting onto 
the star and there is more material in the Keple- 
rian disc which does not project onto the stellar 
disc (surface) for thicker discs, thereby filling the 
central absorption. This is not true for very thin 
discs as it essentially causes the disc and its line 
profile to diminish and both emission and central 
absorption approach the continuum level. Con- 
trary to the case of the wedge-shaped disc (Bu- 
daj & Richards 2004), we observe that in a slab 
disc geometry the blue/red line profile wings of 
the thick discs (a > 3i?0) do not change much 
and emission peaks get slightly closer. 

The middle section of Figure 1 displays the 
effect of changing the outer radius of the disc 
from 6 to 14 Rq. The emission increases as the 
emission volume increases. This parameter has a 
strong effect on the position of the emission peaks 
and width of the central depression. The peaks 
get closer and the central depression narrower for 
larger radius as more matter is involved at lower 
Keplerian velocities. The depth of the central ab- 
sorption varies very little since the amount of the 
matter projected onto the stellar surface does not 
vary as much as in the previous case. This absorp- 
tion first deepens as the overall emission becomes 
more pronounced and then for Rout > 12i?0 weak- 
ens as the fraction of the gas projecting on the 
stellar disc decreases. The blue and red wings of 
the line are not affected at all. 

The bottom part of Figure 1 displays the effect 
of increasing the inner radius of the disc from 2 
to 8 Rq. Contrary to the outer radius, the inner 
radius has little effect on the emission from the be- 
ginning when it is a small fraction of the outer ra- 
dius since it does not change the emission volume 
noticeably. The changes are mainly seen in the far 
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Fig. 1. — Top-left: Effect of varying alpha (vertical half- width of the disc) from 0.5 to 2 Rq; top-right: Effect 
of varying vertical half- width of the disc from 2 to 5 Rq; middle- left: Effect of varying the outer radius of the 
disc from 6 to 9 Rq: middle-right: Effect of varying the outer radius of the disc from 9 to 14 Rq; bottom- left: 
Effect of varying the inner radius of the disc from 2 to 5 Rq; bottom-right: Effect of varying the inner radius 
of the disc from 5 to 8 Rq. 
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Fig. 2. — Top-left: Effect of varying the radius of the star from 1.5 to 3.0 Rq; top-right: Effect of changing the 
inclination of the disc from 90° to 30°; bottom- left: Effect of the geometrical shape of the disc; bottom-right: 
Effect of increasing the microturbulence in the disc from to 60 kms~^. 
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wings of the profile since the higher Kcplcrian ve- 
locities are involved. Only when it approaches the 
outer radius, the overall emission starts to decline 
and the line weakens as the emission volume de- 
creases, and the shape of the profile acquires the 
new U-type shape with a central hole. This shape 
indicates that the central depression is created by 
two different mechanisms (see below). 

In Figure 2 (top left frame), we inspect what 
happens if we vary the stellar radius from 1.5 to 3.0 
Rq . Stellar radius certainly changes the amount of 
matter projected onto the stellar surface but it also 
raises the continumn level and dilutes the emission 
from the disc. This is the main reason why the 
emission weakens for bigger stars. At the same 
time, the central absorption gets wider because 
the matter projected onto the stellar disc spans a 
broader radial velocity interval. When varying the 
stellar radius, it may happen that the star overlaps 
the inner radius of the disc. In such a case, the 
program sets both radii to be equal. 

We also illustrate here the line profiles if the 
disc were viewed from different inclination angles 
(Figure 2, top right frame). Starting from edge-on, 
i = 90°, the two emission peaks have the lowest 
intensity and we also observe the largest separa- 
tion of the emission peaks and the deepest central 
depression. This is because we encounter the high- 
est radial velocities and most of the disc which is 
cooler than the star projects onto the stellar sur- 
face. As the inclination decreases, and we begin to 
see the disc pole-on, the emission peaks increase 
and their separation and central depression van- 
ish, eventually merging into single peak emission 
for i = 0°. One can also see that the central de- 
pression is caused by two different effects. The 
velocity field and geometry of the disc itself can 
produce the U or V type of the depression. Super- 
posed on that is the absorption by the cool mat- 
ter projected onto the hotter stellar surface which 
quickly diminishes if viewed out of the disc plane. 

It is interesting to explore how the shape of the 
profile changes if the geometry of the disc changes. 
The prescription for the velocity field remains the 
same. In Figure 2 (bottom left frame), we compare 
two profiles: one obtained for the slab shaped disc 
with the standard parameters listed above. The 
other disc has the shape of a rotating wedge with 
the wedge angle a = 28°, which produces emis- 
sions of similar strength to the slab disc model. 



All other parameters were the same. It is clear 
that the main difference is in the wings, which are 
more pronounced in the slab model since there is 
more matter closer to the star in that model than 
in the wedge model. 

Apart from the state quantities and velocity 
field, every space point can be assigned a value 
of microturbulence. It is included in the calcula- 
tion of thermal broadening as an additional ther- 
mal motion and can thus be used to model chaotic 
velocity fields on the scale smaller than the mean 
free path of a photon. It turns to be a useful free 
parameter since the mass transfer may not be a 
smooth process and the velocity field of the spi- 
raling gas in the disc may well depart from the 
circular Keplerian orbits and be turbulent, and 
these departures may easily exceed the sound ve- 
locity which is of the order of 10 kms~^. Fig- 
ure 2 (bottom right frame) illustrates the effect 
of increasing the microturbulence in the disc from 
to 60 kms~^. Note that the inner radius Ke- 
plerian disc velocity is about 500 kms"'^. The 
emission peaks are gradually smoothed and are 
broader. The behavior of the central absorption is 
more complicated and interesting. From the be- 
ginning it gets narrower and even deeper because 
of the desaturation effects of the microturbulence 
but then the smoothing effect prevails. The mi- 
croturbulence affects the optical depth along the 
line of sight. While the surface of the star is seen 
at an optical depth of approximately 1.5 (at the 
frequency of the line center) without the micro- 
turbulence, it is seen at the optical depth of 0.7, 
0.4 for the microturbulence of 20 and 40 kms~^, 
respectively. Without the turbulence and assum- 
ing the input parameters mentioned above, the 
line opacity was found about 3 orders of magni- 
tude higher than the continuum opacity in the line 
center. Thomson scattering was the most impor- 
tant continuum opacity source and was a factor of 
300 higher than the hydrogen bound-free contri- 
bution, while H bound-free opacity was a factor 
of 20 higher than the hydrogen free-free opacity, 
and H free-free opacity a factor of 30 higher than 
Rayleigh scattering opacity in the disc. This pat- 
tern, however, will be very sensitive to the state 
quantities, wavelength, and line absorption coeffi- 
cient. 

Figure 3 (top left frame) displays the effect of 
increasing the electron number density of the disc 
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from 1 X lO'^ to 5 X lO^cm"'^. Here, we assume 
that hydrogen is almost Mly ionized, in which case 
the density is hnearly proportional to the electron 
number density. This strongly enhances the emis- 
sion peaks and slightly deepens the central absorp- 
tion. This huge impact on the emission can be 
understood since the equivalent width in the opti- 
cally thin case is proportional to the population of 
the particular level, which in turn is proportional 
to the total HI population. However, the ioniza- 
tion fraction of neutral hydrogen is proportional 
to the electron number density, and if the latter is 
also proportional to the density it follows that the 
total HI population increases with the square of 
the density or the electron number density assum- 
ing the hydrogen abundance and temperature are 
fixed {nHi/riHii ~ UHi/riH ~ uhi/p ~ rig ~ p). 
As a consequence, the equivalent width, EQW, 
should behave also like ~ Uhi ^ n\ which is 
seen in Figure 3. 

Since the emission is so sensitive to the den- 
sity we also studied the effect of different density 
profiles (Figure 3, top right frame). Let the den- 
sity and electron number density be inhomoge- 
neous and have a power law dependence on the 
distance from the star, namely, p ~ r''. This 
was achieved by varying the exponent edendc = r] 
from 0.0 to -2.0 and by normalizing the synthetic 
spectrum to the same emission peak strength. It 
is apparent from the figure that the line profiles 
for lower exponents are broader and have steeper 
wings while those for higher exponents have nar- 
rower and smaller central depression. It can easily 
be explained since the lower exponents emphasize 
the matter close to the star with higher Keplerian 
velocity while the opposite is true for the higher 
exponents. Contrary, to what was found by Budaj 
& Richards (2004) for a wedge-shaped disc with 
the inner radius of 2.5i?0 the position of the emis- 
sion peaks in a slab disc model with = 2.Rq 
is much more sensitive to this exponent. 

Figure 3 (bottom left and bottom right frames) 
displays also the effect of varying the tempera- 
ture of the disc over the interval 5000-20000K. 
The emission is strongest at about 7000K and 
declines towards higher temperatures as the frac- 
tion of neutral hydrogen declines due to ionization. 
The emission also declines towards cooler temper- 
atures as it decreases the population of the lower 
level from which the Ha originates. At the same 



time, the temperature of the disc has a strong ef- 
fect on the depth of the central absorption which 
grows towards the cooler temperatures. For tem- 
peratures T < Tcff, the line source fimction along 
the line of sight hitting the stellar surface steps 
down from the stellar surface to disc and towards 
the observer. For higher temperatures T > T^ff, 
this component of the source function steps up 
reducing the central depression. For both cooler 
and hotter temperatures than 7000K, the line-to- 
continuum emissivity decreases (the continuum is 
provided by the star mainly and is essentially con- 
stant) and the line disappears from the spectrum 
forcing the central depression to merge with con- 
tinuum as well. This is the reason why the central 
depression ceases also below 6000K although the 
line source function steps down. The cooler discs 
if they were more dense could produce even deeper 
central absorption. Nevertheless, one can see that 
discs or regions with temperatures below 5000K or 
above 30000K are not capable of producing strong 
Balmcr line emissions and that most of such emis- 
sions observed in the Algols or CVs must originate 
from the above mentioned temperature regions. 

We would like to emphasize that these calcula- 
tions do not aim to calculate the physical model 
of the disc. The model is assumed to be known. 
The previous calculations illustrate the effect of 
varying just one free parameter at a time while 
keeping all the others fixed. In many cases this 
may not represent the real behavior of the system 
since some quantities may be more or less inter- 
locked such as temperatures and electron number 
densities. 

4.2. Modeling of the Observed Spectra of 
TT Hya 

Now that we understand the influence of var- 
ious parameters on the emerging spectrum, we 
can proceed to fit the complex spectrum of TT 
Hya. We used 23 observed spectra from 1994 of 
Richards & Albright (1999) complemented by 24 
new spectra from 1996/97. A detailed study of 
these spectra will be published in Richards et al. 
(2005) and will include the actual measurements 
of the intensity of the central depression and the 
emission peaks. 

We performed some new and more sophisti- 
cated calculations of synthetic spectra in which 
both stars and a disc were considered. The pri- 
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Table 1: Adopted properties of TT Hya. 
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Designation: u - limb darkening coefficient, a - 
separation of the components, Tp - temperature at 
the rotation pole, /3 - gravity darkening exponent, 
i - inclination, a - vertical half thickness, - 
inner radius. Rout - outer radius, p{Rin) - density, 
TieiRin) - electron number density, rj - exponent of 
the density power law, vtrb - microturbulcncc. Sec- 
ondary was assumed to have synchronous rotation 
and to fill the Roche lobe, VW - Van Hamme & 
Wilson (1993), E - Etzel (1988), L - Lucy (1968), 
BRM - this work. 
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Fig. 4. — Synthetic spectra (solid line) compared 
with the observed spectra (dotted line) of TT Hya 
at different phases. The vertical line indicates the 
laboratory wavelength. The systemic (gamma) ve- 
locity of the system (9.25 kms^-'^) was subtracted 
from the observed spectra. The numbers to the 
right of the spectra indicate the epoch and phase 
of the observations. 
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Fig. 5. — Synthetic (solid line) and observed spec- 
tra (dotted) of TT Hya during total eclipse. The 
numbers to the right of the spectra indicate the 
epoch and phase of the observations. The verti- 
cal line indicates the laboratory wavelength. The 
systemic (gamma) velocity of the system (9.25 
kms~^) was subtracted from the observed spec- 
tra. 



mary was treated as a rotating solid sphere while 
the surface of the secondary was modeled in terms 
of the Roche model. The stellar properties were 
based on the absolute parameters of the system 
obtained by Van Hamme & Wilson (1993) and are 
summarized in Table 1. Since our code accepts the 
polar temperature, Tp, for the Roche model rather 
than the average surface temperature, we adopted 
the value Tp = 4600K which is slightly higher than 
the average temperature T = 4361K obtained by 
Van Hamme & Wilson. Their V-band limb dark- 
ening coefficients were adopted here. The surface 
gravity, log g = 4.23 for the primary and log g = 
2.66 for the secondary were assumed. First, we 
obtained the model atmospheres of both stars by 
interpolating in the Kurucz (1993b) Toff — logy 
grid and assuming solar abundances and a micro- 
turbulence of 2 kms~^. Next, we calculated the 
intrinsic synthetic spectra emerging from these at- 
mospheric models (flux from the unit area of the 
surface) using the code synspec (Hubeny et al. 
1994) modified by Krticka (1998). Here again, so- 
lar abundances and a microturbulence of 2 kms~^ 
were assumed. These synthetic intrinsic spectra 
were then assigned to the primary and secondary 
for the calculation of the complex spectra of both 
stars and the disc using shellspec. Based on the 
lessons learned in the previous section we manipu- 
lated the free parameters of the disc so that the ob- 
served spectra could be reproduced. We assumed 
that the electron number density is equal to the 
hydrogen number density. A good fit to the obser- 
vations was achieved with the realistic parameters 
for the disc. These values are also summarized in 
Table 1. 

Figures 4 and 5 depict both the observed spec- 
tra as well as our synthetic fits for a few represen- 
tative phases at both quadratures and during both 
eclipses to illustrate the most interesting effects. 
The observed spectra display rod and blue emis- 
sion features as well as a central depression. The 
observed far wings of Ha stretching beyond the 
emission peaks in Figure 4 do not seem as dciep as 
in the synthetic spectra. This is not a continuum 
rectification problem. These shallower observed 
wings were previously noted by Peters (1989) and 
interpreted as an emission stretching up to 30A 
from the line center. However, these observed far 
Ha wings in TT Hya differ from single star stan- 
dard spectra or from synthetic ones most probably 
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becaiisc the atmospheric regions of the primary 
are more comphcated. A number of additional ef- 
fects could be expected here like a suction ( "spray 
gun") effect of a disc moving over the atmosphere, 
centrifugal forces, and disc-atmosphere interaction 
with higher turbulence or temperatures. All of 
them tend to reduce the effective gravity of the at- 
mosphere of the primary and the atmosphere-disc 
transition region. Note, that the effective gravity 
is essentially zero in the disc. This could change 
T — T behavior of the atmosphere, mimic lower 
log 3, reduce the pressure broadening, and, conse- 
quently, cause shallower line wings or shell spectra. 
It might be misleading to interpret these shallower 
wings as an emission (Peters 1989). The emission 
peaks and the central depression depend strongly 
on the phase, as best illustrated in Figure 6, for 
all observed and synthetic spectra. A discussion 
of the phase dependence of Ha, simultaneously on 
Figures 4, 5, and 6 follows. 

4-2.1. First Quadrature 

The first quadrature (Figure 4: phase 0.217) 
seems to have the simplest spectrum. Here, we 
see a remarkably good agreement between the 
observed and theoretical profiles. Both emission 
peaks have approximately similar strength. A lit- 
tle bump on the red side of the central depression 
is the intrinsic Ha line of the secondary. 

4-2.2. Secondary Eclipse 

As we proceed towards secondary eclipse, the 
blue emission gets apparently stronger than the 
red emission. This should not happen if the disc is 
symmetric and its source function does not depend 
on the phase. However, this effect is observed also 
in the synthetic spectra (see Figure 6 and Figure 
4: phase 0.475) and is caused by the shadowing of 
the secondary by the disc. Just before secondary 
eclipse the receding parts (red emission) of the disc 
projects onto the surface of the secondary. This 
imposes different boundary conditions for inten- 
sity which is zero in the blue emission but nonzero 
in the red emission. An extra absorption of the 
secondary's light in the receding part of the disc 
causes the red emission to be reduced. This effect 
should be balanced in the middle of the eclipse as 
seen in Figure 6. If this explanation is correct, the 
reverse (stronger red emission) should be observed 
after the middle of the eclipse at phases 0.5-0.6. 



Unfortunately, these phases are not well covered 
by the observations of Richards et al. (2005) but 
Figure 4 of Peters (1989) confirms that this is re- 
ally the case. At the same time during secondary 
eclipse both emission peaks grow slightly. This is 
mainly because the stellar continuum drops dur- 
ing the eclipse and the contribution of the emission 
to the total light can get more pronounced. This 
is also reproduced in the synthetic spectra but the 
observed emission peaks arc still stronger than the 
theoretical ones which indicates that we may have 
underestimated the temperature of the secondary 
at this phase or overestimated the adopted limb 
darkening of the secondary. This might also be 
due to the reflection effect which is not considered 
in our calculations. 

The blue emission extends to higher radial ve- 
locities (Figure 4). This is not predicted by the 
synthetic spectra and our interpretation is that 
this is evidence of the gas stream. At this phase in 
the blue emission we are looking through more in- 
homogeneous regions encompassing the gas stream 
and disc-stream interaction regions which have 
higher than average radial velocity and can be 
rather turbulent exhibiting the strongest depar- 
tures from circular Keplerian velocities. This is 
probably also the reason why the observed central 
depression is displaced slightly to shorter wave- 
lengths. The observed central depression gets very 
deep at these phases, and reaches a minimum at 
about the middle of secondary eclipse. It is deeper 
than in the theoretical spectra. This huge depres- 
sion could be caused by the combination of several 
effects: (1) at this phase the intrinsic Ha lines of 
both stars overlap which deepens the central de- 
pression, (2) the reduction of the emission peaks 
by an additional absorption of light from the sec- 
ondary by the disc might deepen the depression as 
well, (3) as pointed out above, if the surface of the 
secondary were hotter at this phase it would make 
the intrinsic line of the secondary, as well as the 
overall central depression, deeper; (4) the eclipse 
of the secondary lowers the stellar continuum and 
changes both emission and absorption. 

Something unexpected happens after the eclipse 
during phases 0.6-0.7. The blue emission increases 
again. This does not seem to be a temporal ef- 
fect as it can be seen in the earlier observations 
of Peters (1989; Figure 4). One could expect to 
observe the reverse from the opposite line of sight 
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within phases 0.1-0.2 hut it is not the case. Could 
there exist some other circumstellar material apart 
from the disc and the stream in the system that 
contributes to the emission? If this source were 
projected onto the surface of the secondary at 
phases 0.6-0.7 it might have the same effect as the 
disc at phases 0.4-0.5 and might increase the rela- 
tive strength of the blue emission. If this were the 
case, on symmetry grounds, we could expect the 
reverse (red emission stronger) at phases 0.3-0.4. 
There is only a very weak indication in Figure 6 
that this might be the case. Independent data 
of Peters (1989) indicate the same effect that the 
red emission is slightly stronger at phases 0.3- 
0.4. However, if this matter projects onto the 
secondary during 0.6-0.7 and 0.3-0.4 then it would 
have to be outside of the Roche lobe of the pri- 
mary. This may most probably be somewhere 
between the CI and C2 Roche surfaces. 

4-2.3. Second Quadrature 

The emission at the other quadrature (Figure 
4: phase 0.790) looks similar to that at 0.217 and 
is well reproduced in the calculations but now the 
observed central depression is weaker. The central 
absorption is smallest shortly after the quadrature 
at about 0.79-0.86 (Figure 6). Our interpretation 
is that this is evidence of the hot spot or hotter disc 
regions projecting onto the stellar surface of the 
primary which increases the source function along 
the line of sight and reduces the central depression. 
A little bump seen now on the blue side of the 
central depression is the intrinsic secondary Ha 
line. 

4.2.4. Eclipse of the Disc 

The next four phases depicted in Figure 4 
(phases 0.926-0.944) represent the eclipse of the 
approaching part of the disc and one can follow 
in both observed and theoretical spectra how the 
blue emission progressively ceases. The behav- 
ior of both emission peaks is thus explained fairly 
well with our model. Red emission seems to be 
stretching more to the red than predicted but this 
is just the manifestation of the same effect ob- 
served close to the secondary eclipse for the blue 
emission. This time, before the primary eclipse in 
the red emission we observe the gas stream with 
higher receding velocities than expected from cir- 
cular Kepler ian orbits. However, the observed cen- 



tral absorption is again deeper than expected. To- 
tal eclipse which follows is discussed on a separate 
figure (Figure 5) because of the different intensity 
scale. 

The last four spectra in Figure 4 (phases 0.066- 
0.085) cover egress out of the disc eclipse. This 
time, one can observe how the red emission asso- 
ciated with the receding part of the disc emerges. 
Again, the behavior of emission peaks is repro- 
duced fairly well in the calculations, however, the 
observed central depression is much deeper than 
expected from our model and poses the most se- 
rious problem. To some extent it may be caused 
by the similar effect invoked for the explanation 
of the depression near secondary eclipse i.e., (a) 
overlap in radial velocities of stellar components 
(which operates only outside of the total eclipse 
this time). We also encounter (b) the eclipse of 
the low velocity emission rcigions that spill into the 
central absorption which deepens the absorption. 
This is apparently not enough. There are several 
possible explanations of this additional absorption 
just before and after the primary eclipse. It can 
be ascribed (1) to the main body of the disc and 
assumed that the cooler disc regions project onto 
the hot surface of the primary slightly before but, 
mainly, after the total eclipse. We reject this idea 
since such cooler regions in the main body of the 
disc should manifest in asymmetric emission peaks 
at the quadratures which is not observed. An- 
other alternative is (2) to ascribe these deep Ha 
cores to the inner and more dense regions of the 
disc as suggested by Peters (1989). This might 
have no significant effect on the emission peaks 
but there is another problem with this interpre- 
tation and it is the strong phase dependence of 
this feature. There is no obvious reason why this 
absorption if originating from inner disc regions 
should be phase dependent. One could assume 
that an inner disc region close to a hot spot or 
a direct impact of the gas stream is hotter, but 
the matter would quickly cool towards the equi- 
librium value on time-scales much shorten' than a 
Keplerian period of the inner disc. Moreover, we 
observe that an enhanced absorption in the cores 
of Ha are observed only shortly before and after 
the primary eclipse and quickly change and dis- 
appear within 0.1 orbital period between phases 
0.85-0.95 and 0.05-0.15. Such an enhanced ab- 
sorption is observed at the same phases in other 
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Algols e.g., Per (Richards 1993) that docs not 
have a permanent accretion disc. This leads us to 
suggest a third possibility (3) that this effect could 
be associated with sonic cooler matter in the vicin- 
ity of the secondary star or L2 point. One could 
speculate that it may have the form of a slightly 
asymmetric envelope or a "tail" accompanying the 
secondary during its orbital motion or perhaps cir- 
cumstellar matter populating the region between 
the CI and C2 critical Roche surfaces. In either 
situation, it is plausible that such matter could be 
present between CI and C2, and could get there 
easily since the secondary already fills the CI sur- 
face and is a magnetically active star. The disc 
itself reaches almost up to the Roche lobe of the 
primary and there is no reason to believe why it 
should have an abrupt end just before the Roche 
surface and not spill beyond the surface. The lead- 
ing region (in the sense of the orbital movement) 
of such a secondary's envelope should be slightly 
hotter or less dense than the rear or tail of the en- 
velope since the excess absorption is deeper after 
primary eclipse. 

There is still more to understand. Note, that 
just before primary eclipse starts (0.90-0.95) the 
red emission is enhanced while the blue one is be- 
ing eclipsed. The opposite, even more pronounced 
effect is observed after the eclipse (0.05-0.10) when 
the red emission emerges out of eclipse as the blue 
emission gets fainter. It is not clear whether these 
features are permanent or transient or why the 
non-eclipsed emission should be affected during 
the eclipse of the other emission. Possible expla- 
nations include (1) that the disc itself contributes 
significantly to the total continuum and this con- 
tribution is reduced during the disc eclipse or (2) 
the primary is being shadowed by the extra mat- 
ter between CI and C2 (as suggested above). We 
made a test calculation for phase 0.0 with and 
without the disc and found that the contribution 
of the uneclipsed disc to the whole light in the 
continuum at Ha at this phase was only about 
5 X 10^^. Thus, the second option would be more 
plausible. If the effect is real and the explanation 
correct one could expect a drop in the R-band pho- 
tometry light curve at the corresponding phases by 
about 0.05 mag. 



4.2.5. Total Eclipse 

Finally, Figure 5 exhibits the spectra during to- 
tal eclipse. This is an important verification of the 
disc model since at this phase the main source of 
light, the primary star as well as the inner disc re- 
gions, are obscured. The spectral lines of the sec- 
ondary star can be seen more clearly and they indi- 
cate that the temperature of the secondary is prob- 
ably slightly higher than the value of 4361K found 
by Van Hamme & Wilson (1993). We adopted a 
temperature of 4600K at the rotation pole of the 
secondary which is reduced by about 200K on the 
back side of the secondary due to gravity darken- 
ing. These phases also demonstrate that a slight 
decrease in radial density profile of the disc is nec- 
essary to get better agreement with the observa- 
tions. This radial density profile is described with 
the power law behavior p = Pin{r/rin)^ in our 
model and the exponent r] = —1 seems to fit the 
data best. Both emission peaks agree fairly well 
with the synthetic spectra. It is interesting that 
the red emission looks slightly wider than the blue 
one. This is the manifestation of the same effect 
described above at the phases just before primary 
eclipse or of the reversed effect seen close to sec- 
ondary eclipse when the blue emission was wider 
because of the enhanced departures from a Ke- 
plerian velocity field due to the gas stream and 
stream-disc interaction. However, at total eclipse 
the contribution from the remote part of the disc 
close to the Roche lobe and out of the disc plane 
are mainly seen. These areas may exhibit the 
strongest departures from our simple prescription 
of the Keplerian velocity (single star model) which 
might explain why the synthetic central depression 
is wider. 

The calculated central depression during pri- 
mary eclipse is slightly stronger than expected. 
This can be due to several reasons: Firstly, the 
red emission which is wider overlaps the absorp- 
tion and also skews it toward the blue. Secondly, 
the synthetic spectra were calculated for the tem- 
perature of 4600K which may be attained at the 
rotational pole but the actual temperatures we ob- 
serve at the back side of the secondary are lower 
by about 200K and the intrinsic Balmer lines of 
the secondary would be noticeably weaker at such 
temperatures. We would like to note at this point 
that the skew of the central depression to the blue 
at these phases (close to the primary and sec- 
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ondary eclipse) can be easily misinterpreted as ev- 
idence of mass outflow. At the other phases, like 
quadratures, departures from the circular Keple- 
rian orbit may occur causing a shift of the central 
absorption to either side. That is why we feel that 
such blueshifted central depression is not proof of 
the mass outflow from the disc and that this con- 
clusion of Peters (1989) should be revisited. 

4-2.6. Reliability of the free parameters of the 

disc 

Since there are so many free parameters of the 
disc we cannot claim that our fit is a unique solu- 
tion. Nevertheless, as shown in the previous sec- 
tion, the parameters affect the spectrum in specific 
ways that we now understand more clearly (see 
Figures 1-3). Consequently, such calculations can 
be used to put an independent constraint on many 
of these parameters. The calculations are not very 
sensitive to the inner radius of the disc but they 
seem to indicate that it reaches the surface of the 
star, moreover, the density of the disc seems to in- 
crease towards the star. This is expressed by the 
negative exponent rj. However, the exponent is 
best constrained from the emission strength dur- 
ing the total eclipse which is also sensitive to the 
temperature of the secondary and if the latter is 
wrong a spurious value of rj could be obtained. If 
we had adopted the cooler temperature of the sec- 
ondary derived by Van Hamme & Wilson (1993), 
the effect would be to lower the continuum and 
make the emission peaks relatively stronger and 
the exponent rj would have to be even more neg- 
ative. The outer radius can be set much more 
reliably since it affects strongly the separation of 
the emission peaks. The Roche lobe radius of the 
primary is about 1Q.6Rq so our value is realistic 
and is about 95% of the Roche lobe. This is in per- 
fect agreement with the value suggested by Peters 
(1989) who estimated that the disc spans 95% of 
the Roche lobe. The exact agreement is a coinci- 
dence since in practice the disc probably does not 
have an abrupt end but rather dissolves smoothly. 

The vertical dimension of the disc is significant 
and is of the order of the dimension of the primary 
star. It agrees with the estimates of Peters (1989) 
and Plavec (1988) since it is in the region where 
both their estimates overlap. The temperature of 
the disc which we obtained should be viewed as a 
characteristic temperature of the Ha emitting re- 



gions. In reality, there may be strong local depar- 
tures from such a homogeneous model. Regions 
with temperatures between 5000K to 30000K will 
contribute significantly to the Balmer line emis- 
sion but both cooler and hotter regions may man- 
ifest at different wavelengths or spectral lines and 
an image of the disc can change depending on the 
spectral feature studied. The temperature was in- 
ferred mainly from the central depression which 
is strongly variable and there may be other ef- 
fects which are beyond the scope of the paper and 
may turn out to be important like NLTE effects, a 
better treatment of line absorption coefficient and 
scattering at optically thick wavelengths. 

The densities are quite well determined since 
the overall emission is very sensitive to this param- 
eter. Nevertheless, in reality, again strong local de- 
partures can be expected from such a simple den- 
sity profile. The densities we obtained throughout 
the disc decline from rig = 2 x 10^°cm~^ at the in- 
ner radius to rie = 4x lO^cm"^ at the outer radius. 
They are in very good agreement with the estimate 
of Peters (1989; 10^°cm~^). The assumed micro- 
turbulence is very uncertain. In the case of an 
optically thin disc, it should be viewed as depar- 
tures from the Keplerian velocity field rather than 
a local turbulence, and we can claim only that 
these departures will not be higher than about 80 
kms^^. However, these departures may be higher 
in the gas stream. Consequently, we can claim 
that the single star disc model with a Keplerian 
velocity field (as given by Eq. 43) can reproduce 
the shape of the line fairly well. The regions of the 
disc out of the orbital plane and close to the Roche 
lobe are most vulnerable to departures from this 
model. 

The Voigt profile is generally not a good ap- 
proximation for hydrogen lines, however, it can 
be justified in the case of emission lines originat- 
ing from low density regions where the line pro- 
file is determined mainly by the 3D velocity field. 
Finally, we would like to point out that the par- 
allax measured by Hipparcos corresponds to the 
distance of about 154 ± 25pc which is not in very 
good agreement with the previous estimates of Et- 
zel (1988) and Plavec (1988) who both obtained 
distances of about 194 pc. If the new distance is 
confirmed by future parallax measurements, then 
either the geometrical scale (semi-major axis) and 
stellar radii are smaller than those obtained by Et- 
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zcl(1988), Plavcc (1988) or Van Hammc & Wilson 
(1993) or the temperatures of the stars, mainly 
that of the primary, might have to be reduced 
slightly. 

5. Summary and Conclusions 

Wc have developed a new computer code called 
SHELLSPEC. This is to our knowledge the only 
tool which can produce the synthetic light curves 
and synthetic spectra of interacting binaries tak- 
ing into account optically thin moving circumstel- 
lar matter. The code was applied to the well- 
known eclipsing Algol-type binary system TT Hya 
which has an accretion disc. Synthetic spectra of 
the system for all phases were calculated by tak- 
ing into account spherical primary, a Roche lobe 
filling secondary, as well as a disc. For the first 
time, it has been verified by calculations of syn- 
thetic spectra that the double peaked emission in 
TT Hya, and thus also in other long period Algol- 
type eclipsing binaries, originates from circumstel- 
lar matter in the form of an accretion disc with a 
Keplerian velocity field. 

The effect on the spectrum of various free pa- 
rameters was studied. The temperature and incli- 
nation of the disc have the strongest effect on the 
depth of the central depression while the outer ra- 
dius of the disc, the radial density profile, and the 
inclination affect mainly the position and separa- 
tion of the emission peaks. The overall strength 
of the emission is regulated mainly by the den- 
sity and temperature. The central depression was 
found to be created by two different effects: the 
geometry and the velocity field of the disc itself 
and by the cool matter projected onto the hotter 
stellar surface. 

Realistic estimates of the densities, tempera- 
ture, geometry, and dynamics of the disc were 
obtained as a result of these calculations. How- 
ever, there are differences between the observed 
and synthetic spectra which assume a simple cir- 
cular Keplerian disc. These can be attributed 
mainly to the gas stream or stream-disc interac- 
tion. The presence of an additional source of cir- 
cumstellar matter in the system between the CI 
and C2 Roche surfaces is suggested to account for 
the deep absorption in the Ha cores detected near 
primary eclipse and various other observed fea- 
tures. 
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